Distinction between sequential and direct three-body decays 
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We discuss the three-body decay mechanisms of many-body resonances. Sequential decays proceed 
via two-body configurations after emission of the third particle. In direct decay all three particles 
leave simultaneously their interaction regions. The intermediate paths within the interaction regions 
are not observables and only accessible through models. The momentum distributions carry, apart 
from polarization, all possible information about decay modes and resonance structure. In this 
context we discuss detailed results for the decay of the Be(5/2~) resonance. 

PACS numbers: 21.45.-v, 31.15.xj, 25.70.Ef 



Introduction. The distinction between different decay 
mechanisms of a given resonance is a general problem in 
quantum mechanics [l|, l3i Isf • The properties of the frag- 
ments emerging after the decay carry all possible infor- 
mation. A number of accurate and kinematically com- 
plete measurements appear these years in the literature 
for three-body decays [J, S [g]. The results are often 
used to interpret the decay as either sequential via an in- 
termediate configuration or direct to the continuum [7|. 
Also theoretical models are used to classify the decay 
mechanisms [8|, |9| but the results of interpretations from 
different groups of comparable data are not always in 
agreement with each other [3, 0, [l3|. This is unfortu- 
nate since a proper classification of the different decay 
modes is used both in analysis of data and in subsequent 
applications. 

An important question is then to which extent are 
these disagreements real or only related to different mod- 
els and choices of basis. A correct understanding is of 
crucial importance for several reasons, e.g. 

(i) The analyses of data are based on an i?-matrix for- 
malism which intrinsically assumes sequential decay via 
different intermediate configurations. No other proce- 
dure to analyze the data exists. Improvements are in 
addition at best included to account for symmetries and 
the interaction between the first emitted particle and the 
remaining subsystem. 

(ii) The intermediate "paths" leading from initial to 
final states are in quantum mechanics not observables. 
The measured final state information only relates to the 
initial state through the filter of intermediate configura- 
tions. Initial state information and classifications into 
decay modes are then necessarily the results of interpre- 
tations or model computations. Understanding the un- 
derlying physics picture is then necessary for predictions. 

(iii) The interpretations are used e.g. to derive the 
reaction rates for the inverse process of absorption in 
astrophysical environments. Reliable applications must 
use the correct rates which may depend on the adopted 
interpretation of the decay data. 

The purpose of the present letter is to investigate the 



distinct features of different decay modes. We illustrate 
by comparing computed and experimental results and 
conflicting interpretations for decay of ^Be(5/2^). 

Procedure. To describe three-body decay of a many- 
body resonance we must use a three-body model. Effects 
of many-body character at short-distances are mocked 
up by use of a structureless short-range phcnomcnolog- 
ical three-body potential which also is fine-tuning the 
resonance energy. The three particles have masses m^ 
and coordinates ri, i = 1,2,3. To compute the reso- 
nance wave functions we use the Faddeev equations and 
the hyperspherical expansion method ll| combined with 
complex scaling [l^l ■ The crucial coordinate is the hyper- 
radius p defined as a mass averaged distance, i.e. mp^ — 
J2i<j fniruj {ri — TjY /M, where M = mi+ rrik + rrij and 
m is an arbitrary normalization mass chosen to be the 
nucleon mass. The remaining five intrinsic coordinates 
are dimensionless angles. We apply the transformation 
p -^ pexp{i9) for a given scaling angle 9 and compute the 
adiabatic potentials. When 6 > 02 ^ 0.5arctan(r2/2S2) 
where E2 is the energy and r2 is the width, one of the 
potentials describes asymptotically this two-body reso- 
nance with the third particle far away. 

To illustrate we compare measured distributions [7| 
and computations for the 5/2 "-resonance of ^Be (a-l-a-l- 
neutron). The two-body interactions reproduce the low- 
energy a— a and n—a scattering data [i^|lJ| . The three- 
body interaction is used to reproduce the measured reso- 
nance energy above the three-particle breakup threshold 
of 0.856 MeV-zO.4 keV (excitation energy 2.429 MeV). 

Adiabatic potentials. The real parts of the coupled 
set of adiabatic potentials are shown in fig. [TJ The po- 
tential n = 1 has a pocket at small distance holding 
the resonance at the measured energy. All the other 
potentials are repulsive for all distances. The poten- 
tial n = 7 becomes the lowest at around p « 12 fm 
and approaches the ^Bc(0+) resonance energy (dashed 
line) with a corresponding structure combined with the 
neutron far away. The two-body resonance energies of 
^Be(2"'") and ^He(l/2~) are above the first potential for 
almost all distances while the energy of ^He(3/2~) only 



is above the first potential when p is larger than 30 fm. 
The angles 62 for these three resonances are all larger 
than e = 0.1. 
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FIG. 1: Tlie real parts for 9 = 0.1 of the lowest poten- 
tials, including the three-body potential, for the 0.856 MeV 
®Be(5/2~)-resonance (horizontal full line) as function of p. 

The population of the adiabatic components varies 
rapidly at small distances where the coupled adiabatic 
potentials move relative to each other. The couplings, 
when levels get close, are decisive for the choice between 
maintaining the structure or following the energetically 
most favorable path. In the present case about 96% of 
the wave function at large distance maintains the struc- 
ture related to the pocket at small distance. About 3% 
is located in the component related to the ground state 
structure of ^Be. The relative populations become es- 
sentially constant at large distance although the *Be(0+) 
description is less accurate. 

The partial wave decomposition of the dominating com- 
ponent is shown in fig.[21 The angular momentum combi- 
nations at large distance are almost entirely the ^Bc(2+) 
structure of 2h between the two a-particles coupled to 
\h between their center-of-mass and the neutron (fig. [21 
left). The same wave function also has the '^He-structure 
oi Ih between the neutron and one a-particle coupled to 
2h between their center-of-mass and the other a-particle 
(fig. H right). 

Energy distributions. The computed neutron energy 
distribution (left in fig. [3]) almost reproduces the exper- 
imental curve, although lies below at small energy and 
above at high energy. However the data is obviously an 
overestimate at small energy, and due to the normaliza- 
tion this appears as an underestimate at high energy. 
The computed small peak at almost maximum neutron 
energy arises from the potential, corresponding to de- 
cay via the ground state of ^Be. Essentially all other 
parts of the broad and smooth distribution arise from 
the n = 1 potential which in no way corresponds to any 
two-body structures. In fact the distribution solely deter- 
mined from phase space for the largest angular momen- 
tum component in fig. [2] only should be shifted a little 
to match the computed results. The angular distribution 
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FIG. 2: The probabilities \C\ for different angular momen- 
tum combinations {lx,jx,ly,jy) in the two Jacobi coordi- 
nates for the dominating (96%) adiabatic component for the 
^Be(5/2~ )-resonance. 



in fig. [3] between the direction of the neutron and the 
relative momentum of the two a-particles also compares 
favorably with the measurement. 
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FIG. 3: The measured 7\ and computed normalized distribu- 
tions of the neutron energy (left) in units of the maximum, 
and the angle 9i (right) between the relative momenta of the 
Qf-particles and the neutron momentum. Phase space is ab- 
breviated ph.s. 

The traditional analysis in terms of tails of broad 
higher-lying resonances was applied in the recent exper- 
imental investigation f7i| where the *Be(2+) path was 
claimed to account for about 86% of the decay. This 
is in agreement with the theoretical computation in [i3| 
but in contrast to [3| where ^Be(2"'") was ignored and the 
^He-channel was found to contribute a little on top of a 
"democratic" (~ 94%) (phase-space determined) decay. 
In Q it is also claimed that the decay via ^He(p3/2) can- 
not explain the data as well as the ^Be(2+) channel. This 
conclusion is based on analysis of the angular distribution 
while the energy distribution is equally well explained via 
both channels. This difference can be attributed to the 
neglect of ^iie{pi/2) which with ^He(p3/2) is needed to 
be equivalent to ®Be(2"'"), see fig. [21 



Spatial structure The decaying resonance wave func- 
tion depends on the relative coordinates of the three par- 
ticles. The variation describes the dynamical evolution 
of the substructures from before the barrier, via interme- 
diate configurations to large distances. The structure of 
the resonance wave function can be seen in fig. |4] where 
we compare the probability for finding given values of 
the ratio of the distance between two particles and the 
distance from their center-of-mass to the third particle. 

The peak in the first coordinate system (^Be) stabilizes 
at around 0.75 for p larger than about 60 fm. This means 
that when the system emerges from underneath the bar- 
rier at p = 30 fm, see fig. [H the most probable distance 
between the two a-particles is about 18 fm which directly 
shows that the spatial *Be-structure is not present. In 
contrast the peak in the inset moves towards smaller val- 
ues as p increases. This peak arises entirely from the low- 
est potential at large p, see fig.[ll with the two a-particles 
in a relative 0^-structure and a distance of about 3.3 fm 
independent of p. This reflects that ^Be(0+) is popu- 
lated and the a-particles moves apart as for a decaying 
two-body resonance. 

For the small p-values, the other coordinate system 
(^He) reveals a small peak at small values connected to a 
rather broad distribution for higher ratios. The distance 
between the neutron and one a-particle is about 2.8 fm 
for this small peak while the distance to the other a- 
particle is about 27 fm. Combined with the partial wave 
decomposition in fig. [5] this means a small probability for 
a structure resembling ^He. However, as p increases this 
peak disappears before the couplings between different 
potentials have vanished. The distribution stabilizes as 
the broad peak around 1.2 with two a-particles each with 
comparable relative distances to the neutron. 

Sequential and direct decays. The observables are re- 
lated to the asymptotic large-distance results. We clas- 
sify the possible decay modes as direct to the three- 
body continuum or sequential via two-body configura- 
tions. The different decay modes can be found from the 
rotated wave function if the scaling angle is larger than 
the angles of the relevant two-body resonances. 

For direct decays all three particles simultaneously 
leave the interaction region. This corresponds to con- 
verged distributions at the largest accurately computed 
p- values. Further increase of p increases all distances pro- 
portional to p, and the points in fig. [4] can be translated 
into relative distances between the particles. Spatial 
and momentum distributions are asymptotically identi- 
cal since the particles at large distances move along their 
respective coordinates relative to the center-of-mass. 

The decay is sequential if the conserved two-body 
quantum numbers match an intermediate two-body 
structure and the corresponding distance in the complex 
scaled wave function remains small and constant until 
all energy exchanges, or couplings, with the third parti- 
cle can be neglected. Each sequential decay corresponds 
to asymptotic population of the adiabatic component de- 
scribing that structure. These decays result in distribu- 
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FIG. 4; The evolution with hyperradius for two sets of Ja- 
cob! coordinates of the probability distributions for the ra- 
tios between distances or momenta of two particles and their 
center-of-mass and the third particle. The inset shows the 
contribution from n = 7 for *Be. 



tions which continue to change with increasing p like in 
the inset of fig.[H In at least one of the three possible Ja- 
cobi coordinates the distribution is moving towards small 
values of the corresponding distance ratio. Then two par- 
ticles remain close in the intermediate structure while 
the increase of the average size is achieved entirely by 
the third particle moving away. These components must 
be removed from the rotated wave function and treated 
separately as two consecutive two-body decays with pre- 
cisely known kinematics. They are decoupled from each 
other and from the direct decay components. 

One two-body configuration could have attraction 
enough to favor an intermediate resonance structure to 
(exceedingly) large average distances when one particle 
is far outside normal interaction ranges. This could hap- 
pen when the total energy is smaller than the energy of 
the two-body resonance and the effective barrier in fig. [1] 
becomes very wide. The system remains in the classically 
forbidden region to large distances. The signature is sta- 
ble population at intermediate p of one adiabatic com- 
ponent. At larger p the population necessarily decreases 
due to energy conservation corresponding to decay di- 
rectly into the three-body continuum. The three-body 
decay width can be extremely small for such peculiar 
structures which have been baptized as virtual sequential 
decay jl5|. This is decay through the tail of an energy 
forbidden low-lying two-body resonance. The treatment 
can be either to continue the computation to much larger 
p or handle separately as for sequential decay. 

These decay modes are in fact different intermediate 
paths and as such not observables. In pronounced cases 
they leave distinctly different signatures in computed dis- 
tributions as shown in fig. 31 but otherwise a distinction 
is not always meaningful and mode separation becomes 
artificial or model dependent. 



Scaling angle dependence. The presented results are 
all for 6 — 0.1. Increasing 6 to 0.14 and 0.24 beyond the 
^Be(2+) and ^He(p3/2) resonances, respectively, do not 
change the momentum distributions or the fractions of 
sequential decay. However, the accuracy decreases with 
increasing 9 which therefore is kept as small as possible. 

Geometric visualization. One pertinent question is 
whether it is possible to distinguish the different decay 
modes by use of the experimental information without 
theoretical models. Let us plot the measured energy dis- 
tributions P{ei) of all three particles as functions of e^ 
defined as the particle energies divided by the respec- 
tive maximum values. We change P{£i) dci into the dis- 
tribution G{x/y) d{x/y) given as function of the vari- 
able x/y in fig. [4] Since x/y ~ rnijk\/{l — ei)/ei, this 
defines G = 2ei^(l - ei)ej P{ei)/mtjk, where mijk = 
{rrij + mk)y^mJ{mjmkM). 

The resulting plot of G{x/y) is designed to resemble 
the direct decays or the converged part of fig. ID However, 
it contains perhaps also sequential decay contributions 
which may exhibit characteristic features with a peak 
at the corresponding energy and a width equal to the 
sum of the widths of the intermediate two-body state and 
the decaying resonance. For virtual sequential decay the 
peak in the distribution arises from phase space combined 
with the tail of the intermediate resonance. 

To visualize the measured decays geometrically, plots 
of G in fig. [5] could be constructed. One could look for 
features of two-body kinematics corresponding to sequen- 
tial decays. Then it is crucial to correlate consistently the 
results in different Jacobi coordinates. If these features 
cannot be identified it is probably because an extraction 
is ambiguous and the structures can be described in sev- 
eral basis sets or by different models. Then the positions 
of the remaining peaks can be interpreted as length ratios 
in triangular configurations. 

When two particles are identical only two plots are pos- 
sible as for ^Be. The two identical particles go together in 
an intermediate configuration through ^Be(0''') and their 
energy distributions follow from two-body kinematics. 
Everything else is direct decay, and the G— distributions 
in fig. |4] correspond to a triangular decay geometry with 



side ratios d„ 



1.5: 1.3: 1. 



For three identical particles, like ^^C-decay into three 
a-particles, only one such plot exists [9(]. Decays of 
the two lowest 0"'"-resonances are preferentially sequential 
through the ground state of ®Be resulting in a peak cor- 
responding to emission of the first particle and the kine- 
matically related broad structure at lower energy. The 
direct decay of the l+-state gives three peaks and corre- 
sponding triangular side ratios of 2.1 : 1.6 : 1. A directly 
decaying linear chain would produce a distribution with 
peaks at 2/3 and at a very large value. 

Conclusions. The measured energy distributions of 
three-body decaying resonances refiect the structure re- 
lated to the decay mechanism. Specific plots are well 
suited for extraction of the corresponding geometry. 
Sharp peaks most likely correspond to sequential decays 
via energetically allowed two-body configurations. The 
energy forbidden virtual sequential decays leave a signa- 
ture of the corresponding two-body angular momentum 
and parity. Direct decays only carry information of con- 
served quantum numbers of the decaying state. Direct 
and sequential decays can be distinguished if energy or 
angular momentum signatures are different. 

In general, decays via high-lying or broad resonances 
are very unlikely because the corresponding components 
would be depopulated already at small hyperradii where 
the coupling is strong. Sequential decays via energy al- 
lowed low-lying and narrow resonances are possible but 
not unavoidable. Virtual (energy forbidden) sequential 
decays are possible when the two-body resonance energy 
and width both are small and the effective barrier there- 
fore very thick. Several descriptions in terms of different 
basis sets may be equally efficient. 

For application on the three-body decay of ^Be(5/2^) 
we find a small fraction emerging via *Be(0'^) but no de- 
cays via the *Be(2"'") or ^He(p) structures. In this case 
reaction rates based on derived branching ratios via se- 
quential channels are not meaningful or misleading. 
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